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SUMMARY 

2 . SSSSSm £ „ 

eigenvalues and the relationship between line pressure modes and .tom rate mode 
established. A volume at the end of eirchbranch is emp^ J^? Uows my 
combination of boundary conditions, from open to closed, to be used. 

The Jacobi iterative method is used to compute und^ped nanual ^ 

associated pressure/flow modes. Several numerical 

SSEBs for the Helium Supply System of the Space Shuttle Qrbiter Mam 
Propulsion System. 

It should be noted that the method presented hereto can be applied to any one-dimensional 
acoustic system involving an arbitrary number of branches. 

INTRODUCTION 

Often to the analysis of dynamic responses of piped fluid networks, aptrito^-qiuok 
EKaoouT inode shapes and frequencies of the 
mior to the inidadon of more detailed diagnostic team* or 

die fu ndam ental and higher order response mode frequencies o wi»h devices 

SlhSovJfor rapid assessment of modes which may cwmle dynamu^^de™ 8 
^ ^n?3check valves. This provide, valuable dug^^orm ^whe n 

troubleshooting dynamic problems with these types of -5?flS2^v measure and 

and flow mode shapes can provide guidance on positioning of high A^P^^ P^ 8 * 
temdon during testing- Su?h information can be used to infer «MJ^* 
SSSSndSw Stations in regions of the fluid system where 
SJde due to various practical limitations typically encountered on operational system . 

• a 1 _ 


During the course of diagnostic studies of several dynamic p^omenawtorcg^ 
S vK^pellantW systems tmd rocket engines ofthe Space Shuttle*)*, 


iators, 

•iter 
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spacecraft over the last four years, the authors developed a systematic fluid element 
approach to the analysis of related piped fluid networks. These modeling procedures were 
incorporated into a FORTRAN computer program called ACLMODES. 

This paper presents derivations of basic building block equations used in the program, 
illustrates numerical accuracy of the computer code on several problems with known 
closed-form solutions and illustrates how the program was used to analyze several dynamic 
phenomena associated with piped fluid networks of the Orbiter spacecraft 


GOVERNING EQUATIONS 

The basic equations employed in this paper arc the familiar pneumatic/hydraulic 
transmission line equations which govern one-dimensional transient flow. These equa tio ns 
for an unbranched acoustic line are first re-cast in matrix form. Then, this matrix 
formulation is generalized for systems of branched acoustic lines, referred to as "fluid 
networks". 


Matrix Form of Equations for Unbranched Acoustic Lines 

The ordinary differential equations governing one-dimensional flow of an ideal gas, in 
terms of volumetric flow, are (see Figure 1 for notation) 

«i<3i - ^ - Pi+i - Rfi'Qi'Qi O 

q i>l ■ Qj.i • Qj : i-l.N (2) 

where: 

Ij = inertance of the i^ 1 fluid element, 

Cj = capacitance of the i* fluid element, 

0j = mass flow into the i+1 element, 

Pj = pressure at the center of the i* fluid element 
Rfl = resistance, 

N = total number of fluid elements used to model a line segment 

For a uniform line modeled with equal-length elements, the inertance, capacitance, and 
flow resistance are the same for all elements and are given by: 



(3) 

c- ^ 

(4) 

ypRT 


(5) 
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where: 


L = fluid element length 
A = flow area 

Y * polytropic process exponent 

T = temperature 
R = gas constant 

p * density 

f « friction factor (pipe flow) 

D - line internal diameter 
i s equivalent length for minor losses 


a=gaaiS^ g8» !; 

equation of conservation of mass forisentropic flow of an ideal gas. equntion. m 

in fll. [21> &nd r 31 ♦ 


Is 



Figure 1 . Typical discretization of a line segment 


The sets of Eqs. (1) and (2) may be written in matrix form in the special case of R, * 0, that 
is, for the undamped system. The matrix equations are 

A Q + B P = Fj < 6 > 

CP • B T Q = Fj C7) 

where the superscript denotes the transpose of the matrix, and: 


A = 




I 


N 


( 8 ) 
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( 9 ) 


-110 0-00 
0-1 1 0 - 0 0 
0 0-1 1 -oo 


0 0 0 0 -1 1 
0 0 0 0 - 0-1 


( 10 ) 


(B n - -1, B ii+1 - 1, otherwise Bjj - 0) 


Q = < 

r Q ii 

°2 

► ; 

fM 
P 2 > 


A 





( 12 ) 


Differential Equations for Acoustic lines with Branches 

Consider a branched acoustic system such as that shown in Figure 2. The end volumes 
Vj, V 2 and V 3 are used in the formulation for generalizing die boundary conditions. It 
should be noted (see Appendix A) that V * 0 represents a dosed end while V* «* is an 
open end. 
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For simplicity, draping is neglected in this section. 1 Sl^Ete. “ 

^Sn^uTsome"! nodifictians « the ends (elements 

1, K, and N) and element j, where the branch connects to the main line. The first-order 
equations for these special elements are: 


C,Pi = 

Qo-Qi 

(13) 

CJPj - 

Qj-i - Qj - Qk 

(14) 

C k P k = 

Qki 

(15) 

p 

p- 

II 

Qn-i 

(16) 

^ - 

Pj ■ P|C+1 

(17) 


with Cj = 


-Xi_ The B matrix in this case (one branch) has the following structure. 
yRT' 


Bii 


B 

B 


kj 

i i+1 


-1; i = 1, N; i 
-1 

1 ; i - 1. N 


(18) 


Note that B is an N x (N + 1) rectangular matrix. 
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An example on the structure of B: 


j 

k 

N 


2 

4 

7 


(7x8) 


B = 


-1 1 0 0 0 0 0 0 

0-1 1 0 0 0 0 0 

0 0-1 1 0 0 0 0 

0-1 0 0 1 0 0 0 

0 0 0 0-1 1 0 0 

0 0 0 0 0-1 1 0 

0 0 0 0 0 0-1 1 


(19) 


The extension to acoustic lines with multiple branches, or fluid networks, is straight 
forward. Obviously, the structure of B depends on the numbering system used. 


Linearized Form of Governing Equations 
For small pressure/flow oscillations, it can be shown that 

A q + D q + E q q = gj (20) 

C p + H p + E pP = ^ (21) 

where A and C are given by Eqs. (8) and (9) respectively. The diagonal damping matrix D 
is defined by 



where Pj = 2R fi Q® is the linear damping coefficient, Q° bring die mean (steady) flow rate. 

Note that the vectors q, p, fj and f 2 are defined according to Eqs. (11) and (12), with Q, 
P, Fj and F 2 replaced by q, p, fj and f 2 respectively. The matrices H, E q , E p , g t and g 2 
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( 22 ) 


in Eqs. (20) and (21) are given by 

H = B T A 1 D B 1 ' 1 C 

E q * B C 1 B T 
E p « B T A 1 B 

g x = \ - B C * 1 f 2 

g 2 = \ + + DB T * l f 2 ) 


(23) 

(24) 

(25) 

(26) 


ACOUSTIC MODES IN FLUID NETWORKS 
^ . nuid ««wo* - .over^d by 


A q + E Q q = 0 
C p + E p p = 0 


(27) 

(28) 


where A, C, E Q and E p are 


defined by Eqs. (8), (9), (23) and (24) respectively. 


The eigenvalue problem associated with Eq. (27) is 

AX = X q E q X 


(29) 


where Xq = \ and X is a flow eigenvector or flow mode. 

The eigenvalue problem associated with Eq. (28) is 

C Y = XpEpY 

where Xp = \ and Y is a pressure eigenvector or pressure mode. 
(Dp 

It can be easily shown that 

(Dq = ©p 


(30) 


(31) 


and 


Y = C 1 B T X 


(32) 
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Computer Program ACLMODES 

A computer program, referred to as ACLMODES, was developed which computes the 
natural frequencies and associated flow rate/pressure modes for an acoustic line network 
The program can accommodate any number of branches with any combination of boundary 
conditions (ranging from closed to open at each end). The input to die program is relatively 
simple due to its capability of generating acoustic elements with identical properties. 

The program employs the Jacobi Iterative Method to solve the eigenvalue problem defined 
by either Eq. (29) or Eq. (30). Line pressure modes are computed either directly from Eq. 
(30) or using Eq. (32). 


NUMERICAL RESULTS AND DISCUSSION 

Three examples demonstrating capabilities of die ACLMODES program are shown in the 
following paragraphs. The first example is a comparison of a simple system's frequencies 
predicted by ACLMODES and the closed-form solution shown in Appendix A. The other 
examples are actual applications of ACLMODES on Space Shuttle fluid line systems. 


Numerical Test Case 

Example 1 is a test case consisting of a 100 inch long pipe of 0.5 inch LD. filled with 
helium and a volume on both ends. Three different combinations of end volumes, shown 
in Table 1, were used. Note that all volumes are in cubic inches. Case A represents an 
open-closed boundary conditions and Case B a closed-closed. 



mmuzizm 
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Table 1 . Volume sizes for ACLMODES test case. 


ACLMODES was used to determine the first three natural frequencies of each of the three 
cases. Each case was repeated with four different element lengths to evaluate solution 
accuracy versus number of line elements employed. The closed-form solution shown in 
Appendix A was then used to calculated die frequencies of the three cases. The 
ACLMODES results are shown in Table 2 together with die closed-form solution. 
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^HT.T.W faWB 

h risv tt'mm 1 

CASE A 

■ 

100.30 

298.23 

488.85 


1‘HTf 

‘ - ] 1 

iw.w 

301.00 

501.64 

CASEB 


mw 

394.73 

579.92 

200.45 

399.66 

596.40 

200.62 

401.04 

601.06 

2UU.04 

401.24 

601.73 

20U.OJ 

401.30 

601.95 

CASEC 

— I — 
2 
3 

■ 

145.09 

322.36 

511.27 

145.10 

322.89 

513.81 

mm 

145.11 

322.98 

514.28 


Tabic 2. Acoustic Frequencies (Hz) of Cases A, B and C 


•n* perarn difference between the values ACLMODES predicted end dun of the closed 
form solution arc shown in Tabic 3. 



| 


BB 


mokI 

■HH| 




CASE A 

H 

0.10 

0.92 

2.55 

0.02 

0.23 

0.64 

0.00 

0.01 

0.10 

u.uu 

0.01 

0.02 

CASEB 


— 031 

1.64 

3.66 

030 

0.41 

0.92 

0.01 

0.06 

0.15 

OllWJ 

0.01 

0.04 

CASEC 

1 

O.OfM 

0.78 

2.33 

OUT * 

0.19 

0.59 

1 

0.00 

0.01 

0.02 


Tabic 3. Percent Error of Cases A, B and C 
Compared to Closed-Form Solution 


These numerical results show excellent agreement with the ctosed-fam 
ATextJSed? dra U improvement in accuracy of the numencal sotonon as *e number of 
to dSK*s. A general tule of thumb for an acceptable tine element length 

required to obtain accurate numerical results is (oL/c < 0.5, where fi> “ the . • 

circular frequency of the mode sought in rad/sec, L is the line element length (inch) ami c 
the speed of sound in the fluid (inch/sec). Figure 1 shows percent error versus oiVc for all 
values in Table 3. 
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Figure 1. Percent error versus cuL/c 

It can be seen that when the condition coL/c < 0.5 is not met die percent error is greater 
than 1. 


Test Stand Line Dynamics 

Stability testing of the Primary Reaction Control System (PRCS) thruster at the NASA 
White Sands Test Facility (WSTF) required that the test stand have similar line dynamics to 
that of the Space Shuttle Qrbiter. This is because the PRCS thruster is a pressure-fed 
engine so the pressure recovery or waterhammer of the line governs the start-up transients. 
Therefore, a simple line having similar waterhammer characteristics to that of the aft PRCS 
feed system, which is a fairly complicated system with many branches and twelve primary 
thrusters (see Figure 2), was desired. 



Thrusters 


Figure 2. Schematic of the Space Shuttle's PRCS aft fuel supply system. 
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TVe oiimd idea was to use amUar line diameeis to thal of tte 

S^ffrom the supply tank to the thmsters as the test standi™ fcngt^^ver.lt was 
not obvious that this would yield the same dynamics as the vehicle so 
configuration was constructed. 

«E TtSS tntjle jUg***^ ^SiSSiSSSS^ 

S vehide feed 

system. 




Figure 4. First pressure mode of the Space Shuttle s PRCS aft fuel supply system. 
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Space Shuttle Main Propulsion System Helium Supply System 


The regulators in the helium supply system for Space Shuttle Main Propulsion System 
(MPS) were experiencing oscillations. These oscillations were seen both on test stands as 
well as on the vehicle. It was believed that the source of the oscillations was the regulators 
coupling with the downstream line acoustics. 

Each engine has its own helium supply system. A helium supply system consists of high 
pressure supply tanks, tubing leading up to two panels in parallel, and lines from the panels 
rejoining and continuing on to die engine. A panel consist of a regulator and relief valve as 
well as several solenoids and check valves. Only the lines downstream of die regulators 
were of interest, so they were all that was modeled (see Figure 5). 
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. Valve B 


F^elief ^ I 

Sensor B O Regulator B 


Panel B 


V 



Engine 1 


Figure 5. MPS Engine 1 helium schematic of lines downstream of regulators. 


Figure 6 through 8 show the first, fourth and eighteen pressure modes p re di c ted by 
ACLMODES for the Engine 1 helium supply system. The fourth mode is shown again in 
Figure 9 as a flow mode. 
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Figure 6. First pressure mode of the Engine 1 helium 
lines downstream of regulators. 



Figure 7. Fourth pressure mode of the Engine 1 helium 
lines downstream of regulators. 



Figure 8. Eighteenth pressure mode of the Engine 1 helium 
lines downstream of regulators. 
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Figure 9. Fourth flow mode of the Engine 1 helium 
lines downstream of regulators. 


Not all regulators oscillated on the vehicle and those that did, oscillated at different 
frequencies depending on flow demand and number of regulators in use. One mode of 
oscillation was around 1 15 to 120 hertz, with the regulators oscillating out of with 
each other. This mode was predicted by ACLMODES and can be seen in Figure 7 and 
Figure 9. 

The pressure mode shapes were also used to determine if the pressure oscillations being 
measured by a transducer were representative of the oscillations at the regulators. This was 
done by examining the modes with frequencies near die frequency of interest and 
determining if the pressure amplitude at the transducer was being attenuated or amplified 
compared to that of the regulator. 


CONCLUSIONS 

The method presented herein has proven to be a very useful and accurate tool for 
determining dynamic characteristics of complex fluid networks, such as pressure recovery 
and oscillatory behavior. When implemented in a computer code and coupled with a 
plotting routine, this technique can graphically show vital information about the behavior of 
a fluid system impossible to obtain with hand calcnlatinng 
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APPENDIX A 

Modes of a Straight Acoustic Line with General End Conditions 
The classical wave equation for a straight tube in terms of volumetric flow rate is 

(A.1) 

* ax 2 at 2 

The boundary conditions for the system shown below are derived from the continuity 
equation and the definition of fluid capacitance. 



The continuity equation is 



where 

P = 

fluid density 


c * 

acoustic velocity 


Q - 

volumetric flow rate 


p = 

pressure 


A = 

flow area 

But 

ap 
at * 

C (Qin ' Qou.) 


(A.2) 


(A.3) 


where C is the fluid capacitance which is given by 


Gas: 


V _ V _ V 
TfP TRT pc* 


Liquid: 



(B = Bulk Modulus) 
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Substitution of (A.3) in (A.2) yields 


■ - A e«vo»,> 


or, 

* - v(Qi«-Qoui) 


(A.4) 


For end 1, x = 0, 


Qin * ®' Qout B Q I x=0 


Thus, 


SL * V Q| x-0 

Similarly, for end 2, y=L, 

Qin 5 Q I x=L • Qout = 0 


(A.5) 


from which 


ag 

dx 


x=L 



(A.6) 


The general solution of Eq. (A.l) is 

Q(x,t) = T(t) (D,sin ~ x + D 2 cos ® x) (A.7) 


where 


T(t) = BjSincot + B 2 cosox (A.8) 

Substitution of conditions (A.5) and (A.6) into Eq. (A.7) leads to the following frequency 
equation 


£2 tan fl 



'ojn 2 
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< 2-3 


(A.9) 


CO V, "2 

where Q = 7 L is a nondimensional frequency, o. x and Oj « 


Sprtift 1 rjises 
1. V 2 -> 00 (open end) 

1 

Q tan Q = 

a i 


2. V 2 -» 0 (closed end) 

tan Q = - £1 
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